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Abstract 

We calculate long-distance contributions to the amplitudes A{K^ — > 27r, /) in- 
duced by the gluon and the electroweak penguin operators Q% and Qs; respectively. 
We use the l/Nc expansion within the effective chiral lagrangian for pseudoscalar 
mesons. In addition, we adopt a modified prescription for the identification of meson 
momenta in the chiral loop corrections in order to achieve a consistent matching to 
the short-distance part. Our approach leads to an explicit classification of the loop 
diagrams into non-factorizable and factorizable, the scale dependence of the latter 
being absorbed in the low-energy coefficients of the effective theory. Along these lines 
we calculate the one- loop corrections to the 0{p^) term in the chiral expansion of 
both operators. In the numerical results, we obtain moderate corrections to B^^"^^ 

(3/2) 

and a substantial reduction of -Bo . 
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1 Introduction 



In this article we study long-distance contributions to the K ^ tttt decay amplitudes using 
the 1/Nc expansion {Nc being the number of colors) within the framework of the chiral 
effective lagrangian for pseudoscalar mesons. 

The calculation of chiral loop effects motivated by the 1/Nc expansion was introduced 
in Ref. to investigate the A/ = 1/2 selection rule. These articles considered loop 
corrections to the current x current operators Qi and Q2- The gluon penguin operator Qq 
was included at the tree level, consistent with the l/N^ expansion since the short- distance 
(Wilson) coefficient is subleading in N^. Following the same lines of thought the authors 
of Ref. performed a detailed analysis of the ratio s' /e, which measures the direct CP 
violation in K ^ tttt decays. They included the matrix elements of Qq and Qg at the tree 
level in the 1/Nc expansion, arguing that their quadratic dependence on the running mass 
m<j cancels, in the large- A^c limit, the evolution of the coefficient functions in the absence 
of chiral loops. 

In contrast with the A/ = 1/2 rule, which is governed by Qi and Q2, e' /e is dominated 
by the density x density operators Qq and Qs- Therefore it is important to investigate the 
1/Nc corrections to the matrix elements of the last two operators. In particular, it must be 
examined whether the 1/Nc corrections significantly affect the large cancellation between 
the gluon and the electroweak penguin contributions obtained at the tree level in Ref. 0. 
In Ref. the analysis of e'/e was extended by incorporating in part chiral loops for the 
density X density operators, i.e., 1/Nc corrections to the matrix elements of Qq and Qs- The 
final result was an enhancement of {Qq)i=o and a decrease for {Q%)i=2^ which introduces a 
smaller cancellation between these two operators. As a consequence, the authors found a 
large positive value for e'/e [@. 

In this paper we present a new analysis of the hadronic matrix elements of the gluon 
and the electroweak penguin operators in which we include, in the isospin limit, the first 
order corrections in the twofold expansion in powers of external momenta, p, and the 
ratio l/Nci i.e., we present a complete investigation of the matrix elements up to the or- 
ders and /Nc^ One improvement concerns the matching of short- and long-distance 
contributions to the decay amplitudes, by adopting a modified identification of virtual mo- 
menta in the integrals of the chiral loops. To be explicit, we consider the two densities 
in density X density operators to be connected to each other through the exchange of an 
effective color singlet boson, and identify its momentum with the loop integration variable. 

"'^A comprehensive analysis of chiral loop corrections to the 0{p^) matrix elements will be presented 
elsewhere. 
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The effect of this procedure is to modify the loop integrals, which introduces noticeable 
effects in the final results. More important it provides an unambiguous matching of the 
1/Nc expansion in terms of mesons to the QCD expansion in terms of quarks and glu- 
ons. The approach followed here leads to an explicit classification of the diagrams into 
factorizable and non-factorizable. Factorizable loop diagrams refer to the strong sector 
of the theory and give corrections whose scale dependence is absorbed in the renormal- 
ization of the chiral effective lagrangian. The non-factorizable loop diagrams have to be 
matched to the Wilson coefficients and should cancel scale dependences which arise from 
the short- distance expansion. 

The disentanglement of factorizable and non-factorizable contributions is especially 
important for the calculation of the 0{p^/Nc) matrix elements of Qq. although the 
term vanishes for Qq , the non-factorizable loop corrections to this term remain and have to 
be matched to the short- distance part of the amplitudes. These 0{p^/Nc) non-factorizable 
corrections must be considered at the same level, in the twofold expansion, as the 0{p'^) 
tree contribution and have not previously been calculated. The same procedure is followed 
for investigating the matrix elements of Qs- As a final result, we present the numerical 
values for the matrix elements {Q6)o and {Qs)2 to orders and p^/Nc- 

Another improvement is the enlargement of the lagrangian, that is to say, we use the 
complete chiral effective lagrangian up to O(p^). Finally, we include effects of the singlet 
?7o, which is necessary for the investigation of isospin breaking terms. The latter generate 
the matrix element (Q6)2 which is important for e'/e Isospin violating terms will be 
studied in the future. For consistency, and to introduce the general lines of thought, we 
include here the r/o also for the computation of the matrix elements {Q6)o and {Qs)2 in the 
isospin limit, where its effect is expected to be small. 

This paper includes several improvements which are necessary for a complete calcula- 
tion to orders p'^ and p^/Nc, as was defined above. It is still necessary to include these 
improvements for the isospin violating terms, but this will not affect the results for {Qq)o 
and {Qs)2 presented here. Furthermore, we can contemplate still higher order corrections 
which, at present, are beyond the scope of this analysis. 

The paper is organized as follows. In Section ^ we review the general framework of 
the effective low-energy calculation. In Section ^ we discuss the matching of short- and 
long-distance contributions to the decay amplitudes. Then, in Section ^ we investigate the 
factorizable l/N^. corrections to the hadronic matrix elements of Qq and Qs, where we show 
explicitly that the scale dependence resulting from the chiral loop corrections is absorbed 
in the renormalization of the bare couplings, the mesonic wave functions and masses. This 
we do on the particle level, as well as, on the level of the operator evolution for which 
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we apply the background field method. In Section ^ we calculate the non-factorizable 
loop corrections to the hadronic matrix elements and the corresponding non-factorizable 
evolution of the density x density operators. In Section^ we give the numerical values for the 
matrix elements and the parameters Bj^^'^^ and B^^'^\ The latter quantify the deviation of 
the matrix elements from those obtained in the vacuum saturation approximation. Finally, 
we summarize and compare our results with those of the existing analyses. 



2 General Framework 

Within the standard model the calculation of the K ^ ttti decay amplitudes is based on 
the effective low-energy hamiltonian for AS* = 1 transitions 

G ^ 

7if^^=^ = -^enJ]Q(/x)Q,(/i) (/i<m,), (1) 
^ 1=1 

Ciifl) = Ziifi) + TViifl) , T = -it/ in , = V*^Vga , (2) 

where the Wilson coefficient functions Cj(/i) of the local four-fermion operators QiifJ^) are 
obtained by means of the renormalization group equation. They were computed in an 
extensive next-to- leading logarithm analysis by two groups 0, |^. Long-distance contribu- 
tions to the isospin amplitudes Aj are contained in the hadronic matrix elements of the 
bosonized operators, 

(Q,(/i)),^(7r7r, /|g,(^)|ir°), (3) 
which are related to the tt+tt^ and ti^tt^ final states through the isospin decomposition 

(Q,(/x))o = -^(2(7r+7r-|g,(/i)|K°) + (7rV|Q,(/x)|ir°)) , (4) 

{Q^i^^))2 = -^((7r+7r-|g,(/i)|ir°)-(7rV|g,(/.)|X°)) . (5) 

Direct CP violation in K —>■ tttt decays is dominated by the gluon and the electroweak 
penguin operators, i.e., by (^6)0 and (^8)2, respectively, where 

= -2 ^ s(l + 75)gg(l -75)d , = -3 ^ 653(1 + 75)^^(1-75)^, (6) 

q=u,d,s q=u,d,s 

and eg = (2/3, —1/3, —1/3). This property follows from the large imaginary parts of their 
coefficient functions. It is the cancellation between the two penguin contributions which 
gives rise to a small value of the ratio e'/e. Consequently, it is important to investigate 



4 



whether the degree of cancellation is affected by corrections to the hadronic matrix elements 
beyond the vacuum saturation approximation (VSA) 0. 

There are several realizations of non-perturbative QCD [0, |T^, |TT], |T^. A recent de- 
velopment is the calculation of K ^ mi from off-shell K ^ tt amplitudes within chiral 
perturbation theory |jl3|. We will perform our analysis using the 1/Nc approach. To this 
end we start from the chiral effective lagrangian for pseudoscalar mesons which involves 
an expansion in momenta where terms up to 0{p'^) are included ||14|| , 



(7) 



with (A) denoting the trace of A and Ai = dia.g{mu, md, mg). f and r are free parameters 
related to the pion decay constant and to the quark condensate, respectively, with 
r = — 2(gg)//^. In obtaining Eq. (|^) we used the general form of the lagrangian and 
omitted terms of 0{p'^) which do not contribute to the K tttt matrix elements of Qq 
and Qs or are subleading in the expansion.| The fields of the complex matrix U are 
identified with the pseudoscalar meson nonet defined in a non-linear representation: 



U 



exp jU , 



where, in terms of the physical states 

V2K- 



n 



n = vr'^A, , 



v/27r+ 



(AaA;,) — 26ab : 



— TT 



brj' 



\ 



(9) 



gar/ y 



and 



cos6' - V2sin6' , ^26 = sin 6* + ^2 cos 61 . 



(10) 



6 is the 1] — 7]' mixing angle. Note that we treat the singlet as a dynamical degree of 
freedom and include in Eq. (0) a term for the strong anomaly proportional to the instanton 
parameter a. This term gives a non-vanishing mass of the ?7o in the chiral limit (m^ = 0) 
reflecting the explicit breaking of the axial U{1) symmetry. We shall keep the singlet term 
throughout the calculation and will discuss its effects in Section 6. 

^In addition, one might note that the contribution of the contact term cx (Ad'^Ai) vanishes in the isospin 
hmit {niu — md)- 
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The bosonic representation of the quark densities is defined in terms of (functional) 
derivatives: 



r—U^ + L^d^U^d^'UU^ + 2rL8U^MU^ + rH2M^) , (11) 

V 4 / ji 



and the right-handed density {D^)ij is obtained by hermitian conjugation. Eq. ([TT|) allows 
us to express the operators and Qg in terms of the meson fields: 

'1 



q 



1 



+ 0{p'), (12) 



+ 0{p'). (13) 



For the operator the {U'^)(iq{U)qs term which is of 0{p^) vanishes at the tree level. This 
property follows from the unitarity of U . However, when investigating off-shell corrections 
it must be included. This important aspect, which was not studied previously, will be 
discussed in detail in the following sections. 

The 1/Nc corrections to the matrix elements {Qi)i are calculated by chiral loop dia- 
grams. The diagrams are ultraviolet divergent and are regularized by a finite cutoff. This 
procedure, which was introduced in Ref. |]l[, is necessary in order to restrict the chiral la- 
grangian to the low-energy domain. Since we truncate the effective theory to pseudoscalar 
mesons, the cutoff has to be taken at or, preferably, below the 0{1 GeV). This limitation 
is a common feature of the various phenomenological approaches, which at present do not 
include higher resonances. 

The loop expansion of the matrix elements is a series in 1//^ ~ l/^o which is in direct 
correspondence with the short-distance expansion in terms of a^/vr ~ l/Ai^^: the large- A^c 
behaviour of SU (Nc) quantum chromodynamics is represented by diagrams which have a 
planar gluon structure. Subleading terms in the 1/Nc expansion are included by means of 
internal fermion loops (suppressed by a factor 1/Nc) or non-planar gluon interactions (sup- 
pressed by [0. These corrections actually generate the multimeson intermediate 
states which constitute the loop diagrams of the effective theory.^ 

A pedagogical introduction to the 1 /Nc expansion in terms of mesonic degrees of freedom can be found 
in Ref. 111. 
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Finally, we note that the meson loop corrections are needed not only for improving 
the matching of the matrix elements to the short- distance coefficient functions but also 
for obtaining the correct infrared structure, which is required to maintain the unitarity 



relations at low energy [17, 18 



3 Matching of Long and Short Distance 

To calculate the amplitudes we follow the lines of Ref. |l|] and identify the ultraviolet cutoff 
of the long-distance terms with the short-distance renormalization scale fi. In carrying out 
this matching we pay special attention to the definition of the momenta inside the loop. 
This question must be addressed because the loop integrals, within the cutoff regularization, 
are not momentum translation invariant. 

In the existing studies of the hadronic matrix elements the color singlet boson connect- 
ing the two densities (or currents) was integrated out from the beginning ||, ||, ||. Thus 
the integration variable was taken to be the momentum of the meson in the loop, and the 
cutoff was the upper limit of its momentum. As there is no corresponding quantity in the 
short- distance part, in this treatment of the integrals there is no clear matching with QCD. 

The ambiguity is removed, for non-factorizable diagrams, by considering the two den- 
sities to be connected to each other through the exchange of the color singlet boson, as 



was already discussed in Ref. ||18[. A consistent matching is then obtained by assigning the 
same momentum to the color singlet boson at long and short distances and by identifying 
this momentum with the loop integration variable. This important feature of the modified 
approach is illustrated in Fig. 1. The momentum of the virtual meson is shifted by the 
external momentum, which affects both the ultraviolet, as well as, the infrared structure of 
the 1/Nc corrections. The same method was used in studies of the K^ — Ks mass difference 
||T9| and the evolution of current x current operators in the chiral limit . 



Obviously, the modified procedure described above is not applicable to the factorizable 
part of the interaction. However, in the next section we will show explicitly that all fac- 
torizable terms quadratic and logarithmic in the cutoff are independent of the momentum 
prescription in the loop. Moreover, they are absorbed in the renormalization of the bare 
low-energy coefficients, as well as the mesonic wave functions and masses. Consequently, 
the factorizable 1/Nc corrections are not to be matched to any short-distance contribu- 
tion, i.e., they refer to the strong sector of the theory. Therefore there is no need for a 
momentum cutoff, and we will calculate the remaining finite corrections using dimensional 
regularization, which constitutes a momentum invariant procedure. 
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4 Factorizable 1/Nc Corrections 



Since factorizable and non-factorizable corrections refer to disconnected sectors of the the- 
ory (strong and weak sectors), we introduce two different scales: Ac is the cutoff for the 
factorizable diagrams and Ac for the non-factorizable. We will refer to them as the fac- 
torizable and the non-factorizable scales, respectively. A similar analysis of chiral loop 



corrections was performed to determine the parameter ^1 



We shall prove in this section, within the cutoff regularization, that the quadratic and 
logarithmic dependence on Ac which arises from the factorizable loop diagrams is absorbed 
in the renormalization of the low-energy lagrangian. Consequently, in the factorizable 
sector the chiral loop corrections do not induce ultraviolet divergent terms, i.e., the only 
remaining ultraviolet structure of the matrix elements is contained in the overall factor 
~ I /ml. This is to be expected as the evolution of m^, which already appears at leading 
A'c, is the inverse of the evolution of a quark density. Therefore, except for the scale of l/m^ 
which exactly cancels the factorizable evolution of the density x density operators at short 
distances, the only scale remaining in the matrix elements is the non-factorizable scale 
Ac. It represents the non-trivial part of the factorization scale in the operator product 
expansion. Since the cutoff Ac disappears through renormalization, the only matching 
between long- and short-distance contributions is obtained by identifying the cutoff scale 
Ac of the non-factorizable diagrams with the QCD renormalization scale /z. 

The proof of the absorption of the factorizable scale Ac will be carried out in the isospin 
limit. This explicit demonstration is instructive for several reasons. First, we verify the 
validity of the general concept in the case of bosonized densities which, contrary to the 
currents, do not obey conservation laws. Second, we check, within the cutoff formalism, 
whether there is a dependence on a given momentum shift {q ^ q±p). Thirdly, including 
the rjo as a dynamical degree of freedom we examine the corresponding modifications in 
the renormalization procedure. Finally, there remain finite terms from the factorizable 
1/A^c corrections which explicitly enter the numerical analysis of the matrix elements. This 
point will be discussed at the end of this section. 



4.1 Calculation of the Matrix Elements 

Due to the unitarity of the matrix field U the tree level expansion of Qq starts at the O(p^). 
Consequently, including only the first order corrections in the twofold expansion in external 
momenta and the ratio 1/A^c, no additional terms arise from the renormalization of the 
wave functions and masses, as well as, the bare decay constant / since these corrections 
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will be of higher order. This statement does not hold for the electroweak operator Qs 
which, for tt+tt^, induces a non- vanishing tree matrix element at the 0{p^). 

The wave function and mass renormalizations can be deduced from the pion and kaon 
self-energies, i.e., from a calculation of the propagators at next-to- leading order in the 

1 /2 

double series expansion. For the wave functions we obtain (defining iir = Zj ttq) 



1 , 8-^5 2 
1 H — —m^ 



P 



+ 



(47r)2/2 3(47r)V^ 



m 



K 



+ 



2ml 



3(4vr)2/^ 



log ( 1 + 4 



mt 



1 + -^ml - 



A 



+ 



\og\l 1 



^ ^ ^m^ + 2ml) + ■ 



(14) 
(15) 



8Z/5 2 

1 + -j^m^ 



+ 



(477)2/2 4(47r)2/2 



A 



ml log 



+ cos^^ m; log ( 1 + ) + sin^^ m;, log [ 1 + 




1 + 



8L,, 



-m 



K 



logA^ 1/^2 2\ 

+ TT^^7:{^rn\ + ml) + 



(47r)2/2 (47r)2/2 6 



2m| log 1 + 



A? 



(16) 



(17) 



where the ellipses denote finite terms we omit here for the analysis of the ultraviolet 
behaviour. One might note that Eqs. (|T^ and (IB) are exact only if the cutoff is associated 
to the virtual meson. However, any momentum shift (g ^ g ± p) modifies only the finite 
corrections (compare Eq. (|69|) of Appendix B). 

In specifying Eq. ([T7|) we applied the octet-singlet squared mass matrix. 



4m\ — ml 



-2V2{ 



m 



K 



ml 



m 



-2^/2{mj^ - ml) 2m^ + ml + 3a 



with a 



m 



li + ral^i — 2m]^ and the corresponding mixing angle p2 

3a 



tan 26* 



2mlo 



2V2 



m\ 



00 



mt 



1 - 



2(m 



mt 



(19) 



The mass renormalization is found to be 

,2 



rm 



2 (L5 - 2Ls) + -a 



/ 



3 (47r)2/2 



m + ms 



m 



K 



r- 



im 



K 



(L5 - 2L8) + -a 



+ • 

log A 



3 (47r)2/2 



+ 



(20) 
(21) 
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where rh = {rriu + md)/2. The ratio of Eqs. ( pO[ ) and (|2T|), to one-loop order, determines 
the difference L5 — 2^8 of the low-energy couplings: 



m 



K 



mt 



rh + Trig 
2m 

m + ms 
2m 



[m 



K 



m 



P 



^{L,-2Ls 



+ ■ 



1 - 



(m 



K 



m 



F2 



(22) 
(23) 



Note that Eq. (^) exhibits no explicit dependence on the scale Ac; i.e., the chiral loop 
corrections of Eqs. ( pO|) and (pTf ) do not contribute to the SU (3) breaking in the masses 
and, consequently, can be absorbed in r. This implies (modulo finite terms) 



-^5 ~" 2Lg 



Finally, / and L5 are obtained from the decay constants of pions and kaons ^ 



1 H —mt — 



1 + 



P 
4L, 
P 



+ 



2{47iyp ' {47ryP2 
3 



-m 



K 



+ T-r^T7TT7 7(5m^ + mj 



2{An)^P (477)2/^4 



+ ■■ 



+ 



(24) 

(25) 
(26) 



Defining the constant through the relation 



-1^ = 1 + 



4L| 
F3 



(m^ - ml) 



(27) 



from Eqs. (pSf ) and 



we find, to one-loop order 
L5 



jr 



3 1okA2 

^ + 



16 {Any 



(2J 



which is in accordance with the result from chiral perturbation theory |T1|. Then, from 
Eq. (g4|) we get 

3 log A^ 



jr 

^8 



32(4.)^+-- 

One might note that the coefficient in front of the logarithm in Eq. ( pQ]) differs from the 
one given in Ref. |]14|. This property follows from the presence of the singlet tiq in the 
calculation. Eqs. (p2|) and ( p3|) define the renormalization conditions because the term 



I/5 — 2I/g plus the constant terms which appear in the ratio of the masses in Eq. (^ 



determine the bare constant L5 — 2Ls. Similarly Eqs. (p5|)-(p7D with the associated finite 
terms determine the coupling constant L5. As we focus in this section on the ultraviolet 
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behaviour we omit the finite contributions. Full expressions relevant for the numerical 
analysis are given in terms of integrals in Appendix A. 



Next we investigate the (bare) tree level of the K'^ — >■ tttt matrix elements, up to 0{p'^) 
in the chiral expansion, as well as, the factorizable 1/Nc corrections to the 0{p^). The 
latter corrections refer to the first term on the right-hand side of Eqs. ([T^) and ([T3|) . Both 
contributions can be calculated from the diagrams depicted in Fig. 2. From the sum of 
these diagrams we obtain 



^(7r°7r°|Q6|^") 

l{Ti^Tl~\Q,\K'') 

z(7r°7r°|g8|K^) 
z(7r+7r-|Q8|ir^) 



0\F 
(0) 



0\F 
(0) 



0\F 
(0) 



4v^ 
"T 

T 



2 / 2 

r ( — m 



2/2 2 - 



-r (m^ — 



^5 + 



3 logAg 
16 (47r)2 

3 logA^ 



16 (47r)2 

3 log Ag 
16 (47r)2 



+ 



+ 



+ ■ 



(30) 
(31) 
(32) 



0\F 
(0) 



4 



1 -T77^(^i + 2^')(^5 -121^8) - 



3A? 



3/ 



(47r)V^ 



+ 47^?^(21m^ + 24m^ + 8a) 



+ ■■ 



(33) 



12 (4vr)2/2 

The structure of Eqs. (|30|)-(^2|) guarantees that the renormalization of L5 renders them 
finite. The situation is more complicated for the matrix element in Eq. (|33D as we will 
remark below. 

If we use Eqs. (|l5l)-(^), including only the first order corrections in the parameter 
expansion, we arrive at the renormalized (factorizable) matrix elements of the Qq and 
operators:"^ 



0\F 
(r) 



4^2 / 2ml 



K 



m + nis 



[m 



K 



(r) 



0\F 

(r) 



4^2 / 2m]^ 



m + TUs 



2V2 ( 



2m 



K 



K 



(34) 
(35) 

(36) 



z{n+n-\Qs\KY 



(r) 



2m 



K 



\m + m^ 



2 r 



3F^ 



im 



K - 11^')^5 



"^Lg does not appear in the matrix elements of Qe because its contributions to the first and second 
diagram of Fig. 2 are canceled by a tadpole (third diagram of Fig. 2). 
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1 fi 

~{m\-2ml)Ll 



+ •••• (37) 



Eqs. (|3^) -(^7|) are unambiguous, as the quadratic and logarithmic terms in Eqs. (p!5|) -(^3l) 
were found to be independent of the momentum prescription inside the loops. 



Note that the factorizable scale Ac is absent in Eqs. (|3^ -(p7D [except for the running of 
l/(m + ms)^ ^ 1/m^]. Residual scale dependences could nevertheless unfold at the orders 
p^/N^ or p^/Nc- The latter would arise, e.g., if we used / rather than F.^ in the O(p^) tree 
level expressions of Eqs. ( p3D and (^) or Eqs. (p4D-(p^). This would be consistent at the 
level of the first order corrections in the twofold series expansion, as the difference concerns 
higher order effects. However, the scale dependence of / (which is mainly quadratic) will 
be absorbed by factorizable loop corrections to the matrix elements at the next order in 
the parameter expansion and has not to be matched to any short-distance contribution. 
Consequently, it is a more adequate choice to use the physical decay constant in the 
expressions under consideration. Instead of F^r the kaon decay constant Fx could be used 
as well. Both choices will be considered in the numerical analysis, which gives a rough 
estimate of higher order corrections. At the same level of accuracy, in the O(p^) terms of 
Eqs. (|3^ -(p7D the prefactor [2m|;/(m + m^)]^ could be replaced by (m^/m)^. However, 
this choice is unsuitable as m suffers from larger uncertainties. 

We note that the coefficients in front of and Lg in the matrix element of Eq. ( pT]) 
are different from those of the bare couplings L5 and Lg in Eq. (|33|) . The change of 
the coefficients comes about as the quantities in Eq. (|33|) are replaced by renormalized 
quantities. In particular, the quadratic term in Ac is absorbed in the renormalization of 
the decay constant / and the mesonic wave functions. Finally, omitting the constant terms 
which refer to the factorizable loop corrections, Eqs. ( |5B| ) and ( P?] ) are combined to obtain 
the isospin-two tree level matrix element of Qs up to 0{p'^) in the chiral expansion: 



271 V 



2 \ 2 



tree ^ j K 

2 I ^ _|_ 



4 o „ o. 16 

Fr 



F^ + —{2ml - ^ml)Ll - — (m^ - 2m^)L^ 



(3J 



The numerical value for this term is given in Table 1. In Ref. ||2^ only the bare matrix 
elements were included in the corresponding tree level analysis of {Q8)2- Consequently, the 
new contribution of Ref. pSl, i.e., the contribution of the Lg coupling, was found with a 



sign opposite to that in Eq. (PB]). This was corrected in Ref. ||2^ in the framework of the 
chiral quark model. 
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4.2 Operator Evolution 

The results of the previous section can also be seen directly at the operator level, in par- 
ticular at the level of the density operator. To demonstrate this we apply the background 
field method as used in Refs. and for current x current operators. This approach is 
powerful as it keeps track of the chiral structure in the loop corrections. It is particularly 
useful to study the ultraviolet behaviour of the theory. 

In order to calculate the evolution of the density operator we decompose the matrix U 
in the classical field U and the quantum fluctuation ^, 

U = exp(^e//) U, e = TAa 

with U satisfying the equation of motion 

Ud^Lf^ - d^tJU^ + rtJM^ - rMU^ = ^(In f/ - In f/^) ■ 1 , U = exp(i7r"Aa//) 



(39) 
(40) 



The lagrangian of Eq. (|^) thus reads 

The corresponding expansion of the meson density around the classical field yields 

The evolution of {DLjij is determined by the diagrams of Fig. 4, and we obtain 

logA^ 



(41) 
(42) 



4r(t/t),.(0) + ^r(f/t),.(0)^ 



12^ '^ ^ ' (47r)2 



H2 + 



3 log A 
16 "(4^ 
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L8 + 



3 log a: 

32 (47r) 



21 



L5 + 



3 log A 



21 



(43) 



16 (47r)2 

The quadratic and logarithmic terms for the wave function and mass renormalizations can 
be calculated from the diagrams of Figs. 6 and 7, i.e., from the off-shell corrections to 
the kinetic and the mass operator, respectively, second and third term of Eq. (^if ). The 
resulting expressions for and mj^ turn out to be identical to those found in the explicit 
calculation of the pion and kaon self-energies, Eqs. (pO]) and (pi]). For the wave functions 
we get 



8L 



P 
8L5 



5 2 



log A 



(47r)2/2 2 -(47r)2/ 



2 f2 ' 



-m 



K 



'(4vr)2f 



—m 



logA^ 



K 



(44) 
(45) 
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Along the same lines and Fk can be calculated, to one-loop order, from the diagrams 
of Fig. 5, and we obtain^ 



■ ^ 4L5 2 3 A, 
1 H — —m^ — 



1 + 



AL 



logA^ 



2(47r)2/^ 4 "(47r)2/2_ 



r 



5 2 



H — m 



logA^ 



2(47r)2/^ ' 4-^^(47r)2/ 



(46) 
(47) 



In accordance with Eqs. (p!5|)-(|33|) both the quadratic and the logarithmic terms of 
Eqs. (^3l)-(0) prove to be independent of the way we define the integration variable inside 
the loops. This is due to the fact that the quadratically divergent integrals resulting from 
the diagrams of Figs. 4-7 [i.e., those of the form d'^q/ {q±p)^ ] do not induce subleading log- 
arithms, that is to say, all quadratic and logarithmic dependence on the scale Ac originates 
from the leading divergence of a given integral. 

Now looking at Eqs. 



-(^Tf) we observe that the ratio 11// and, consequently, the 
matrix field U are not renormalized (i.e., ttq// = vr^/F^ and Kq/J = K^/Fk)- Then, by 
means of Eqs. (^) and (^61) , we can rewrite the density of Eq. (|^) as 



2m\ 



(m + Ills 



F2 

TT 



8L| 
F2 



K 



ml) 



F2 



-m 



K 



(4J 



with X = diag(m^, m^, 2m|; 



In obtaining Eq. 



couplings of Eqs. (ESI) and (E^). In addition, we introduced 



+ 16^4^ + 



we used the renormalized 



(49) 



Note that the renormalized density exhibits no dependence on the scale Ac, except for 
the scale of l/(m + ms). Note also that in Eqs. (^3]) and ( ^51) we did not specify logarithmic 
terms induced at the one- loop order which correspond to the L4, Lg and L7 operators in 
the chiral effective lagrangian of Ref. [0]. An explicit calculation of these terms shows 
that they give no contribution to the K tttt matrix elements of Qq and Qg. 

In conclusion, using a cutoff regularization the factorizable contributions to the and 
Qs operators up to the orders and p^/Nc are given, modulo finite loop corrections, in 
terms of the K ^ tttt matrix elements by Eqs. (^)-(|37|) or in terms of a single density by 
Eq. (^81). Our results exhibit no explicit scale dependence. Moreover, they do not depend 



^The representation of the bosonized current in terms of the background field can be found in Ref. [po|. 
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on the momentum prescription inside the loops. The finite terms, on the other hand, will 
not be absorbed completely in the renormalization of the various parameters. This can be 
seen, e.g., from the fact that the rescattering diagrams of Fig. 2 contain a non- vanishing 
imaginary part. As the renormalized parameters are defined to be real, the latter will 
remain. 

In addition, the real part of the finite corrections carries a dependence on the momentum 
prescription used to define the cutoff. However, we proved that the chiral loop diagrams do 
not induce ultraviolet divergent terms. Therefore we are allowed to calculate the remaining 
finite corrections in dimensional regularization, which is momentum translation invariant 
(i.e., we are allowed to take the limit Ac — > oo). This procedure implies an extrapolation of 
the low-energy effective theory for terms of 0(m^_^/A^; m'^^^/X'^; . . .) up to scales where 
these terms are negligible. This is the usual assumption made in chiral perturbation theory 
for three flavors. 



5 Non-factorizable 1/Nc Corrections 

The non-factorizable l/N^. corrections to the hadronic matrix elements constitute the part 
to be matched to the short-distance Wilson coefficient functions; i.e., the corresponding 
scale Ac has to be identified with the renormalization scale n of QCD. We perform this 
identification, as we argued in Section |, by associating the cutoff to the effective color 
singlet boson. Then, at the 0{p^) in the chiral expansion of the Qq and Qs operators, from 
the diagrams of Fig. 3 we obtain 

|Q,|A-»)- ^ 4 (J^)^ - ™i) + . . . , (51) 



,V.-|Q,|A-°r = + (53) 



Again, for the reason of brevity in Eqs. (|50|)-(p3D we did not specify the finite terms which 
must be included in the numerical analysis (in particular, they provide the reference scale 
for the logarithms). In addition, we replaced m^, m^, and the mixing angle 6 by m^, 
and the instanton parameter a using the octet-singlet mass matrix of Eq. (0). 
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Note that in Eqs. (|50|)-(|53D we used and 2m\/{rh + m^) rather than the bare 

parameters 1// and r. Again the difference represents higher order effects. However, the 
appearance of / or r in Eqs. (^D])-(^) would induce a dependence on the factorizable scale 
Ac, which has no counterpart in the short-distance domain (compare the discussion at the 
end of Section |4.1|) . As for the factorizable contributions the choice of instead of 
would be also appropriate. 

The results presented above are in accordance with the non-factorizable evolution of Qe 
and Qs we obtain in the background field approach by calculating the diagrams of Fig. 8: 



Q^F^|^2^ r.2 I '^'nT'K ^ ' log A 



rh + nis J (47r 



(54) 



+^{d,U^UMUd^'U% + ^(f/tx + x^U)dsSgg + ^a(f/t),,(f/),. 



• (55) 



Only the diagonal evolution of Qe, i.e., the first term on the right-hand side of Eq. (|5^), 
gives a non-zero contribution to the matrix elements of Eqs. (^D|) and (pT]). In particular, 
the mass term which is of the Lg and H2 form vanishes for K ^ tttt decays, as do the 
Lg and H2 contributions at the tree level (see Section In Eq. (|55|) for completeness we 
kept the terms proportional to 6qq which, however, cancel through the summation over the 
flavor index. 

Note that Eqs. (|5^ ) and (^) are given in terms of operators and, consequently, can 
be applied to K ^ Sn decays, too. Note also that our results, Eqs. (p0|)-(|55|), exhibit no 
quadratic dependence on the scale Ac; i.e., up to the first order corrections in the twofold 
series expansion the matching involves only logarithmic terms from both the short- and 
the long-distance evolution of the four-quark operators. This is due to the fact that there 
is no quadratically divergent diagram in Fig. 8 apart from the first one which vanishes for 
the Qq and Qs operators. Moreover, for a general density x density operator there are no 
logarithms which are the subleading logs of quadratically divergent terms. Therefore, all 
the logarithms appearing in Eqs. (|50|)-(|55D are leading divergences, which are independent 
of the momentum prescription. The finite terms calculated along with these logarithms 
depend on the momentum prescription. They are, however, uniquely determined through 
the matching condition with QCD (see Fig. 1). 
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One might note that the statements we made above do not hold for current x current 
operators: the l/N^. corrections to these operators, performed in the first non- vanishing 
order of their chiral expansion, exhibit terms which are quadratic in Ac. Furthermore, 
already these terms were shown to depend on the momentum prescription |P0|| . 

We close this section by giving the long-distance evolution, at the 0{p^), of a general 
density X density operator Q'^^^'^ = —8{Dji)ab{DL)cd- As we showed in Section the 



factorizable 1/Nc corrections do not affect its ultraviolet behaviour. Then, from the non- 
factorizable diagrams of Fig. 8 we find 



/^abcd 



+ - 



(0) 



■2a logA2 



A? 



da jrbc 



2m\ 
171 + rris 



log A 



(4vr) 



{U^X + X^Uf^S^" + 5'^''{xU^ + Ux^f 



(56) 



The corresponding expressions for the non-factorizable loop corrections to operators Q 
and Qs; Eqs. (ISl) and (pH), can be obtained directly from Eq. 



6 Numerical Analysis and Discussion of Results 

To compute the hadronic matrix elements of Qq and Qs we pursued the following strategy. 
First, the non-factorizable contributions were calculated, in the isospin limit, from the 
diagrams of Fig. 3. In this part of the analysis the finite terms, which are systematically 
determined by the momentum prescription of Fig. 1, are completely taken into account. By 
using algebraic relations all integrals resulting from the diagrams of Fig. 3 can be reduced to 
three basic integrals which are given explicitly, in the framework of the cutoff regularization, 
in Appendix B. They were computed up to terms of the order and p^, respectively, that 
is to say, to a relative precision of approximately 10~^. Second, the finite terms arising 
from the factorizable loop diagrams of Fig. 2, as well as, from the renormalization of the 
wave functions, the masses and the low-energy couplings were estimated using dimensional 



regularization, as discussed at the end of Section ^ 

We use the following numerical values for the parameters: 

= (m^o+m^+)/2 = 137.3 MeV, = 92.4 MeV, 

TTiK = (m^o+m^+)/2 = 495.7 MeV, Fk = 113 MeV, 

= 547.5 MeV, 9 = -19°, 

m^/ = 957.8 MeV . 
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Substituting them in Eqs. (|2|) and we compute LI = 2.07x10-^ and = 1.09 x 10^^ 
For the numerical values given above, is close to 2Lg, and we find the 0{p^) tree 
level contribution to {Qs)2 to be small, because the term proportional to m\ in Eq. ( PS]) 
approximately vanishes. This result is different from the statements made in Ref. [^]. The 
full expressions needed for the renormalization of the parameters /, L5 and Lg in terms of 
integrals are presented in Appendix A. In the octet limit the results in the appendix are 
the same as in Refs. |]14| and Finally, we used the ratio nis/rh = 24.4± 1.5 which 



enters in the calculation of Lg. 

The values we obtain for the matrix elements of Qq and Qs are presented in Table 1, 
where we also specify the various contributions coming from the factorizable and the non- 
factorizable diagrams, respectively. In these matrix elements we have extracted the factor 

= [2m|^/(m + m^)]^, whose dependence on the factorizable scale will be canceled ex- 
actly, for a general density x density operator, by the diagonal evolution of the Wilson 
coefficients. Finally, for comparison, we present in Table 1 also the numerical values ob- 
tained by replacing F^^ by in the 0{p^) and 0{p^/Nc) factorizable and non-factorizable 
contributions, that is to say, in the corresponding terms of Eqs. (0), ( p7D [or Eqs. (|60D-(|63D 
of Appendix A], (^)- (|37|) and (^OD-(^), and in the finite terms. The difference gives a 
rough estimate of the higher order corrections. 

In Table 2 we list the corresponding values for the Bi factors, which quantify the 
deviation of the hadronic matrix elements from the VSA results: 

= \{Q,)o/{Q6)r^\ , Bf'^ = \{Qs)2/{Qs)r''\ ■ (57) 

The VSA expressions for the matrix elements were taken from Eqs. (XIX. 16) and (XIX. 24) 
of Ref. m. Numerically, they are |(Q6)o| = 35.2 ■ i?^ MeV and |(g8)2| = 56.6 MeV ■ [/^^ _ 
(0.389 GeV)^]. The second term in the expression for Qg contributes at the 2% level and 
has been neglected in Table 2. 

We discuss next the corrections to the matrix elements {Q6)o and {Q8)2- As already 
mentioned, the operator Qq is special because the 0{p^) tree level matrix element is zero 
due to the unitarity of the matrix U. Nevertheless the one-loop corrections to this matrix 
element must be computed. These corrections are of 0{p^ /N^ and are non-vanishing. 
We have shown in Eqs. ( |50l) and (|5T|) that the explicit calculation of the loops yields a 
cutoff dependence (i.e., a non-trivial scale dependence) from the non-factorizable diagrams 
which has to be matched to the short-distance contribution. In addition, the logarithms 

^Note that our constants L5 and Lg should not be confused with the scale dependent coefficients L5 
and LI in Refs. ITl and 123. 
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Ac = 0.6 GeV 


A, = 0.7 GeV 


A, = 0.8 GeV 


A, = 0.9 GeV 


• \ NF loops 


-35.2 
-68.4 - 37.0i 
29.8 + 37.0i 


-35.2 
-68.4 - 37.0i 
34.6 + 37.0i 


-35.2 
-68.4 - 37.0i 
39.0 + 37.0? 


-35.2 
-68.4 -37.0i 
42.9 + 37.0i 




38.6 
(45.8) 


33.7 
(41.8) 


29.4 

(38.2) 


25.5 
(35.0) 


• / ^ \ NF loops 


56.4 
56.0 -O.lz 
-20.7- n.5i 


56.4 
56.0 -O.li 
-24.8 - 11.5i 


56.4 
56.0 -O.H 
-28.8 - 11.5i 


56.4 
56.0 -O.H 
-32.8 - 11. 5z 




37.2 
(40.2) 


33.2 
(37.0) 


29.5 

(33.8) 


25.9 
(30.7) 



Table 1: Hadronic matrix elements of Qq and Qs (in units of R^- MeV), shown for various 
values of Ac. The numbers in the parentheses are obtained by replacing F^^ by Fx in the 
next-to-leading order expressions. 

of Eqs. ( PH ) and ( P3| ) are needed in order to cancel the scale dependence of various bare 
parameters in the tree level expressions, as was checked explicitly in Section ^. We note 
that in the twofold expansion in and p^/Nc the contribution of the loops over the 
0{p^) matrix element must be treated at the same level as the leading non- vanishing tree 
contribution proportional to L5. This is revealed by the large size of the non-factorizable 
0{p^/Nc) corrections presented in Table 1. It is the sum of both, the factorizable and the 
non-factorizable contributions, which renders the numerical values for (^6)0 close to the 
VSA value. For the imaginary part, which is due to on-shell rescattering effects and does 
not depend on the matching condition with QCD (see Fig. 1), the cancellation is complete. 
This property follows from the {U'^)dq{U)qs structure of the operator. The main effect of 
the loop corrections is to change the dependence of {Qg)o on Ac, from a flat behaviour at 
the tree level to the dependence presented in Tables 1 and 2, which is important for the 
matching. We note that at Ac ~ 700 MeV the value for the matrix element of Qq is very 
close to the VSA result leading to i?6 — 1- 
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Ae = 0.6 GeV 


A, = 0.7 GeV 


A, = 0.8 GeV 


A, = 0.9 GeV 


^(1/2) 


1.10 


0.96 


0.84 


0.72 




(1.30) 


(1.19) 


(1.09) 


(0.99) 


o(3/2) 
-"8 


0.66 


0.59 


0.52 


0.46 




(0.71) 


(0.65) 


(0.60) 


(0.54) 



Table 2: and Bg factors for various values of the cutoff Ag. 



The Qs operator is not chirally suppressed, i.e., its 0{p^) tree level matrix element is 
non-zero. In this article we include the tree level contribution up to O(p^), as well as, 
the loop corrections of 0{p^/Nc), that is to say, corrections to the first term of Eq. ([I3|). 
This is a full leading plus next-to- leading order analysis of the Qs matrix element. The 
one- loop corrections, even though suppressed by a factor l/N^ with respect to the leading 
tree level, are found to be large and negative, leading to the small values for Bg presented 
in Table 2. These large corrections persist in the octet limit [i.e, in the absence of the r/o, 
with a = b = 1 and = (4m|^ — m^)/3]. Therefore, they are not due to the presence of 
the rjo which brings in a small change. One might note that the numbers in Table 2 are 
specified for the central value of rris/m = 24.4 ± 1.5 Including the error of this mass 
ratio changes the Bg parameter by ± 0.06. 

In comparison with {Qe)o, the non-factorizable corrections to (^8)2 are less pronounced, 
as expected from the power counting scheme in and 1/A^c- However, because the fac- 
torizable 0(j5^) and 0{p^/Nc) corrections are small (and negative), the non-factorizable 
terms produce a significant reduction of {Q8)2- The size of the higher order terms indicates 
that the leading- Ai'c calculation or the closely related VSA are not sufficient for the matrix 
elements of the Qs operator.j] 

In view of the large corrections one might question the convergence of the 1/A^c ex- 
pansion. However, there is no strong reason for such doubts because the non-factorizable 
contribution we consider in this paper represents the first term in a new type of a series 
absent in the large- Nc limit. It is reasonable to assume that this leading non-factorizable 
term carries a large fraction of the whole contribution. 

(3/2) 

As a general result, we note that our study indicates a significant reduction of -Bg . By 
comparison the corrections to B^^^"^^ are moderate, i.e., there is no clear tendency for values 

^This has already been observed for the matrix elements of Qi and Q2 which are relevant for the 
A/ = 1/2 selection rule ^. 
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much larger or smaller than one. Our results for (Q6)o and {Qs)2 can still be improved by 
calculating higher order terms in p"^ and 1/Nc, like for instance those of 0{p'^/Nc) which 
will be along the lines of this work. The 0{p'^/Nc) will bring in a quadratic dependence 
on Ac ||29[ and even though suppressed by a factor of p"^ relative to the 0{p^ /Nc) may 
compensate, to a large extent, the scale dependence of the logarithmic terms of this paper. 
Another improvement would be to include the vector mesons which is a new calculation 
beyond the scope of this work. 

It is interesting to compare our results with those of other calculations. Refs. and 
01 investigated 1/Nc corrections to the matrix elements of Qq and Qs- This calculation 
considered the product of the two densities without the color singlet boson and the match- 
ing of short- and long-distance contributions was not explicit as in the present article. The 
0{p^/Nc) contribution to Qq was not included, but terms of 0{p'^/Nc) were included in 
Qe and Qg. In this study the parametrization of the 0{p^) lagrangian was not general 
as only one coupling constant was introduced. The numerical results showed also a ten- 
dency of reducing {Qs)2 substantially, whereas {Qq)o was found to be enhanced compared 
to the VSA result. Calculations in lattice QCD obtain values for Bq close to the VSA, 
sJ^/^^(2GeV) = 1.0 ± 0.2 0, ^ and 0.76(3) (5) |1. Recent values reported for Bg are 
sf/^^(2GeV) = 0.81(3)(3) 0.77(4)(4) g, and 1.03(3) These studies use tree 

level chiral perturbation theory to relate the matrix elements {7T7T\Qi\K) to (irlQilK) which 
are calculated on the lattice. The chiral quark model |^| yields a range for Bq which is 
above the VSA value, -Bg^''^^(0.8 GeV) = 1.2 — 1.9, and predicts a small reduction of the Bs 
factor, i?g^^^^(0.8 GeV) = 0.91 — 0.94. Although the scales used in the lattice calculations 
and the phenomeno logical approaches are different, the various results can be compared as 
the Bq and B^ parameters were shown in QCD to depend only very weakly on the renor- 
malization scale for values above 1 GeV [0 . Finally, in their analysis of e'/ e the authors 
of Ref. considered Bq and Bs as free parameters to be varied around the VSA values 

o{l/2) _ p(3/2) _ 

We note that our result for B^^"^^ is in rough agreement with those of the various 

(3 /2) 

studies quoted above, whereas the value we obtain for B^ lies below the values reported 
previously. It is desirable to investigate whether this substantial reduction of {Qs)2, which 
is due to the non-factorizable l/Nc corrections to the leading term in the chiral expansion 
of Qs, will be affected by higher order corrections. This point is of great phenomenological 
interest because a less effective cancellation between the Qq and operators, in the range 
obtained in the present analysis, will lead to a large value of e' /e in the ball park of 
~ 10~^. This can be confirmed or disproved by the forthcoming experiments at CERN 
(NA48), Fermilab (E832) and Frascati (KLOE). 
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7 Summary 



It was recognized, long ago, that the operators Qq and Qs are of central importance for 
the determination of the CP parameter e'/e. This makes the calculation of their matrix 
elements imperative as the Wilson coefficients are known to a good degree of accuracy. 
We carried out this calculation in the 1/Nc expansion, where we included terms up to 
and 0{p^/Nc). In doing so we introduced several improvements. First we used 
the complete pscudoscalar lagrangian relevant to these orders and included effects of the 
singlet rjQ, which we found to be small. At the same time we paid special attention on 
the definition of the momenta in the chiral loop corrections. To this end, we considered 
the exchange of a bosonic field between the quark densities whose momentum is taken to 
be the same at long and short distances. In this approach we set up the identification 
of the ultraviolet cutoff of the long-distance terms with the QCD renormahzation scale. 
This procedure leads naturally to the classification of the diagrams into factorizable and 
non-factorizable. 

We showed explicitly, to C(p°/iVc), that the factorizable scale of the chiral loop cor- 
rections is absorbed in the renormahzation of the low-energy lagrangian. Thus for the 
factorizable terms the matching of long- and short-distance contributions is between the 
running quark masses and quark densities where the matching is exact, i.e., the scale de- 
pendence drops out. There remain the non-factorizable terms where we showed explicitly 
that the dependence on the cutoff, to the order of our calculation, is only logarithmic. Our 
analysis was carried through using two different techniques. The first one is an explicit 
calculation of the matrix elements at the particle level which involves a large number of 
diagrams. The second is the background field method. It leads to operator relations which 
can be applied also to X — > Stt decays. We verified that both techniques give the same 
results for the quadratic and logarithmic terms. The full finite corrections were calculated 
at the particle level. 

Finally, we determined the numerical values of the matrix elements. We obtained 
moderate corrections to (Q6)o and a large decrease of {Qs)2- Each of these matrix elements 
depends on the renormahzation scale, but it is significant to emphasize that their ratio is 
fairly stable. The numerical results indicate that loop corrections are important and must 
be included. We note that the terms of 0{p^/Nc) we calculated here arc lowest order 
corrections to the well established 0{p'^) chiral lagrangian and do not contain any large 
ambiguity. It remains to be seen whether the results of Tables 1 and 2 will be affected by 
higher order corrections. This point is important because a cancellation between the gluon 
and the electroweak penguins in the range obtained in the present analysis will lead to a 
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large value of e' je ~ 10 
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A Bare Parameters 



In terms of the basic integrals the full expressions needed for the renormalization procedure 
read 



8L 1 / 

SL 1 / 

Zk = 1 + 'JT^'k - (yh [^tt] + 2/1 [rriK] + cos^O h [m^] + sin^6' Ji [m^/ 

8?7i^ 1 / \ 

1 - ^(i^5 - 2L8) + — (^3/iK] - a2/i[m,] - 262/i[m^,]j 



= rm 



m 



m + nis 



K 



r- 



1 



(58) 
(59) 
(60) 



5m 



/ 



f{L,-2Ls 



Ii[m^] (m;(a^ - 46^) - 8m|f(a - 6)6 - m^(a + 26)^ 



+2/i[m^/] 2m^a(a + 26) - m^(a^ - 6^) - mi {a - 6) 



l + ^m2-^(2/iK]+/i[mx] 



1 + -j^^K - (lilmn] + 2/1 [mx] + cos^6' /i[m^] + sin^6' /i[m^' 



(61) 
(62) 
■ (63) 



a, b and 6* are defined in Eqs. (|TOD and (p^Qf), the integral Ii[m] in Eq. (0) of Appendix B. 
Eqs. (p9D-(|6l|) and (^) differ from the corresponding expressions in Ref. ||14| on account 
of the presence of the rjo state. In the octet limit [9 = 0, = (4m|^ — m^)/3] Eqs. (pB[)- 
(|63|) are in agreement with Ref. We note that the rjQ state modifies the logarithmic 

dependence of the Lg coefficient on the renormalization scale, whereas it does not affect 
the corresponding term in L5. 

In the cutoff regularization scheme Eqs. (|58|)-(|63D together with the explicit form of 



the integral Ji given in Eq. (|67|) of Appendix B lead to the formulas listed in Section 
in which the finite terms have been omitted. 

In the numerical analysis of the matrix elements dimensional regularization has been 
used for the factorizable part, as argued at the end of Section |4.2| , and the integral Ji has 
been calculated in the standard way. The full expressions for the (renormalized) parameters 



^The comparison is carried out by omitting the L4 and Lq terms which are subleading in Nc- 
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/, Ls and Lg have been obtained from Eqs. (|60|)-(p3|) by replacing / byF^^ (or Fx) in the 
C(p^) and OijP /Nc) terms, as discussed at the end of Section 



B Basic Integrals 

Using algebraic relations the complex structures of the four-dimensional integration can 
be reduced to three basic components: 

Performing a Wick-rotation to Euclidian space-time the ultraviolet cutoff may be imple- 
mented through the step-function 6'(A^ — g|). A straightforward calculation then yields 

1 

16^ 

In order to determine I2 and Is analytically we shift the variable q by the external mo- 
mentum. Properly taking into account the resulting modification of the upper bound we 
introduce an angular-dependent argument of the step-function. However, we omit the 
explicit angular-integration writing the latter function in terms of a Taylor- series: 

e(M -,S^a)^ .(A? - + . (68, 

The corresponding solution of the integral I2 reads 



Ii[m] 



rn? 



(67) 



/AV /AX [,2 3 , .0. 8, 



The computation of requires a Feynman-parametrization: 

h[mi,m2,p\ = dx J ^^{(g-xp)^- [xV-a;(p^ + m^-m^)+m^]}"^ (70) 
Performing the Wick-rotation and introducing the variable q'^ = Qe — we obtain 

Islmi, 1712, Pe] = -jI^ j ^^'^''' [q'^ + M^{x)f ^ ~ ^'^ ~ '^xWp)e - x^pl] , (71) 
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with 

M^(x) = —x'^pI + x{pI — ml + ml) + m\ . 
For distinct masses mi and m2 Eq. ([7l|) yields 

a2" 



(72) 



/a [mi, ma, p] 



167r2 



arctan 



2 2 

p m2 



2{m\ — mly 



m2 
mi 

1 



1 + log 



m^ — m^ + 

mima 
Al + ml 

2^2 



+ arctan 



+ 



m^ 



m; 



_(A2 + m2)2 



(2A^m^ + mfm2 + ma) + 



ms 



2m? 



log 



— ml + p2 
+ ml 



2 2 

m| — m2 



log 



A2 + mf 
A? 



m? 



Al + ml 



4 2 

p m2 



mt 



m 



2\4 



6(A2 + m|)4 {^^c'fnliml + ml) + 3A^^ml{-ml + Gmjml + 7m\ 



+2A^m2(2m]^m2 + ITm^^mg — ma) + m2{m^ + lOmimg + mg 



ml{ml + ml) f + ml 



[mf — ml) 
where we defined 



log 



A2 + m2 



+ 0(p' 



/ 2 , 2 2\2 ^22 

w = [m^ + ma — p ) — 4m^m2 • 



(73) 



(74) 



As Is starts only logarithmically in the cutoff dependence, in Eq. (|73D we truncated the 
series including only terms up to the order p^. 

In the case of equal masses we perform a power series expansion with respect to the 
parameter 5m^ = mi — ml- Then putting 5m^ to zero we find 



h[mi = m2,p] 



167r2 



p2 



arctan 



p^ 



mt 



A2 + m? 



log 



m; 



A2 + ml 



2 (3Ag + m?)m2 ^ (-20A^ + SA^m? + mf)m2 

^/ A 2 , ~~2V5 ^ P 



60(A2 + m2)5 



+ 0(p^). (75) 



6(A2 + m2)3 

with w being reduced to w = p*^ — Ap'^mf. 

Through analytic continuation, Eqs. (^) and (^) provide complex solutions. These 
appear for a/p^ > mi + m2. In the process under consideration, the latter relation can only 
be satisfied for mi = m2 = m-j^, p = px- Thus our analysis gives the physical threshold 
condition for tt — tt rescattering, the imaginary part of being attributed to the strong 
final state interaction phase. 
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Figure Captions 



Fig. 1 Matching of short- and long-distance contributions. 

Fig. 2 Tree plus factorizable loop diagrams for the K ^ mr matrix elements of Qq and 
Qs] the crossed circles denote the bosonized densities, the black circles the strong 
interaction vertices. The external lines represent all possible configurations of the 
kaon and pion fields. 

Fig. 3 Same as in Fig. 2, now for the non-factorizable loop diagrams. 

Fig. 4 Evolution of the density operator; the black circle, square and triangle denote the 
kinetic, mass and Ua{^) breaking terms in Eq. (^l]), the crossed circle the density of 
Eq. (^21). The lines represent the C, propagators. 

Fig. 5 Evolution of the current operator. The crossed circle here denotes the bosonized 
current. 

Fig. 6 Evolution of the kinetic operator (wave function renormalization) . 
Fig. 7 Evolution of the mass operator (mass renormalization). 

Fig. 8 Non-factorizable loop diagrams for the evolution of a density x density operator. 



30 




Fig. 3 
31 





Fig. 4 



Fig. 5 



Fig. 6 



Fig. 7 




Fig. 8 



